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We propose a theoretical scheme to show the possibility ol generating motional nonlinear coherent 
states and their superposition for an undamped vibrating micromechanical membrane inside an 
optical cavity. The scheme is based on an intensity-dependent coupling of the membrane to the 
radiation pressure field. We show that if the cavity field is initially prepared in a Fock state, 
the motional state of the membrane may evolve from vacuum state to a special type of nonlinear 
coherent states. By examining the nonclassical properties of the generated state of the membrane, 
including the quadrature squeezing and the sub-Poissonian statistics, we find that by varying the 
Lamb-Dicke parameter and the membrane's reflectivity one can effectively control those properties. 
In addition, the scheme offers the possibility of generating various types of the so-called nonlinear 
multicomponent Schrodinger cat sates of the membrane. We also examine the effect of the damping 
of the cavity field on the motional state of the membrane. 

PACS numbers: 42.50.Wk, 42.50.Ct, 42.50.Dv 
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I. INTRODUCTION 

The coupling of mechanical motion of mcsoscopic or 
macroscopic mirrors to the electromagnetic degrees of 
freedom via radiation pressure[U [2] is a promising ap- 
proach for applications ranging from the detection of 
weak forces and small displacements to fundamental 
studies of the transition between the quantum and the 
classical world [5] . The canonical optomechanical system 
consists of an optical cavity where one of the end mir- 
rors is free to move [4j. The radiation pressure acting 
on the movable mirror realizes an optomechanical cou- 
pling between the incident optical modes and the various 
vibrational modes of the mirror. The use of optome- 
chanical coupling has been proposed many years ago for 
the implementation of quantum limited measurements 
of mechanical forces [5], as the interferometric detection 
of gravitational waves |6j or atomic force microscopy [7]. 
Since then, optomechanical systems have generated much 
experimental and theoretical interest [BHT2] . They offer 
the prospect of realizing quantum effects at a macro- 
scopic scale |13j . of supplying quantum sensors for ap- 
plications ranging from single molecule detection |14j to 
gravitational wave interferometry [151 116] , for the quan- 
tum control of atomic, molecular, and optical systems 
j!7j . and for possible new quantum information process- 
ing devices 18J . The optomechanical cooling and trap- 
ping of mirrors has also recently become the subject of 
an intense research effort [5 HTT1 fT51 IT5] as it offers a viable 
means of exploring quantum effects such as superposition 
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and entanglement at a macroscopic scale |13l |2"U] . 

The main technical barrier to reaching the regime in 
which either the mechanical device itself or its readout 
demonstrate quantum behavior, has been difficulty of 
integrating ultrasensitive micromechanical devices with 
high-finesse optical cavities. Recently, it has been de- 
veloped a novel type of optomechanical system, the so- 
called mcmbranc-in-thc middle geometry, for addressing 
this issue in which the mechanical degree of freedom is a 
flexible, partially transparent dielectric membrane placed 
inside a Fabry-Perot cavity with fixed end mirrors [2"lT 
[23] . This has the advantage of not having to combine 
the flexibility needed for the mechanical oscillator with 
the rigidity of a high-finesse cavity mirror. Although the 
membrane is nearly transparent, it couples to the opti- 
cal cavity dispersively. This coupling is strong enough to 
laser-cool a 50-nm-thick dielectric membrane from room 
temperature(294 K) down to 7mK|21j. In addition, the 
dispersive nature of the optomechanical coupling allows 
for realization of a sensitive displacement squared read- 
out of the membrane|21j. Such a readout is a crucial re- 
quirement for measuring quantum jumps in a mechanical 
oscillator. It has been shown theoretically[24] that an op- 
tomechanical system with the membrane-in-the-middle 
geometry can be used to detect phonon Fock states. Fur- 
thermore, for that system it has been predicted [25] two- 
phonon cooling of the mechanical oscillator, squeezing of 
the mechanical oscillator, and squeezing of the optical 
output field. The extension of cavity optomechanics to 
multimcmbrane systems has also been considered [26 . 

In connection with quantum state engineering, op- 
tomechanical systems have also attracted considerable 
attention because of the great possibilities they have to 
produce nonclassical states of both, the mechanical oscil- 
lator and cavity fieldgl E Ql EH [23 UZHSO]- Similar to 
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the case of a light field where the imprint of the quantum 
regime is signaled by generation of nonclassical states, the 
generation of nonclassical states of a mechanical oscilla- 
tor can be a hallmark for quantum control of a macro- 
scopic massive object[3T]. The possibilities to generate 
nonclassical states, and in particular a superposition of 
coherent states of the quantized field are because of the 
Kerr-like Hamiltonian that may be obtained in the op- 
tomechanical systems [23 HI] • Besides, the fact that the 
motion of the mechanical oscillator is quantized allows 
the possibility to generate a large variety of nonclassical 
states of the mirror [2"?]. It is shown that the mirror can 
also be prepared in a Schrodinger cat state with many 
components by a quadrature measurement of the cavity 
field after its interaction with the moving mirror. More- 
over, the possibility of generating even and odd coherent 
states of a moving mirror has been considered [32]. A de- 
tailed analysis of the effects of finite temperature on cre- 
ating and verifying a macroscopic quantum superposition 
in a micro-optomechanical system has been presented in 
[33j . It has been shown that an unambiguous demonstra- 
tion of a quantum superposition requires the mechanical 
resonator to be in or near the ground state. In Ref.|8 
the authors propose a scheme for transferring quantum 
states from the propagating light fields to macroscopic, 
collective vibrational degree of freedom of a massive mir- 
ror by exploiting radiation pressure effects. Some of other 
more recent proposals include the generation of squeezed 
states of a vibrating membrane or a movable mirror via 
reservoir engineering in an optomechanical setup [50]. the 
creation of EPR beams in an optomechanical system [51]. 
and the creation of robust entanglement against increas- 
ing temperature between the optical intracavity mode 
and the mechanical mode of the mirror [55]. 

In recent years, there has been paid much attention 
to the study of a family of generalized coherent states, 
the so-called nonlinear coherent states(NLCSs)[36] be- 
cause of their relevance in nonlinear quantum optics. 
They correspond to nonlinear algebras rather than Lie 
algebras[37]. For a detailed discussion of the math- 
ematical properties of NLCSs we refer the reader to 
Ref.[3"Tl |3"8] . The NLCSs, and their superposition ex- 
hibit nonclassical features like amplitude squeezing, sub- 
Poissonian statistics and self-splitting accompanied by 
pronounced quantum interference effects [39 . These 
states are not mere mathematical objects. In the con- 
text of ion trap, the microscopic NLCSs may appear 
as stationary states of the center-of-mass motion of a 
trapped and laser-driven ion far from the Lamb-Dicke 
regime 36J . The method makes use of the strong nonlin- 
earities inherent in the Jaynes-Cummings Hamiltonian 
for the laser-assisted vibronic interaction [40 . Very re- 
cently, it has been introduced [H] a physical scheme that 
allows one to prepare and control a special kind of these 
states, the so-called sphere-coherent states, as the mo- 
tional dark states of a properly laser-driven trapped-atom 
system. Furthermore, it has been explored 42 the possi- 
bility of generating various families of NLCSs in a lossless 



micromaser within the framework of intensity-dependent 
Jaynes-Cummings model. 



Motivated by the above-mentioned studies on the gen- 
eration of various types of nonclassical states in an op- 
tomechanical structure, in the present paper we deal with 
the question of how the NLCSs can be generated in an 
optomechanical system. The system under consideration 
is an optomechanical cavity with membrane-in-the mid- 
dle geometry that consists of a high finesse cavity with 
two perfectly reflecting fixed end mirrors, and a par- 
tially reflective movable middle mirror (such as a dielec- 
tric membrane) [2"TH2"3l |4l5] . In this type of optomechani- 
cal structure, the radiation pressure and cavity detuning 
are periodic in the membrane displacement. This pe- 
riodicity leads to an intensity-dependent interaction of 
the radiation pressure field with the movable membrane. 
As we shall see, in the absence of any damping mecha- 
nisms, this system can lead to the generation of motional 
nonlinear coherent states and their superposition for the 
membrane. We also study the nonclassical properties 
of generated state and show that how the system pa- 
rameters (in particular, membrane's reflectivity and the 
Lamb-Dicke parameter) can strongly affect the dynam- 
ics of the system under consideration. If those states are 
to be observed, then very good isolation of the system 
from the environment is necessary. Although, by tak- 
ing into account the relevant experimental parameters 
ED HI], just the same considered in this paper, one can 
neglect the membrane's motional damping, but one can- 
not neglect the cavity damping. Therefore, in order to 
make the model under consideration more realistic, we 
also consider the photon leakage from the cavity as a rel- 
evant source of decoherence and examine its influence on 
the motional state of the membrane. 



The paper is organized as follows. In Sec. II we first 
derive an intensity-dependent Hamiltonian describing the 
coupling of the micromcchanical membrane to the radi- 
ation pressure field through j-phonon excitations of the 
vibrational sideband. Then, by restricting our treatment 
to the case of j=l and by using the Feynman disentan- 
glement theorem we derive the time evolution operator 
of the system. In Sec. Ill we consider the generation of 
NLCS of the membrane and investigate the nonclassical 
properties of the generated state. In Sec. IV the prepa- 
ration of the nonlinear multicomponent Schrodinger cat 
states of the membrane is discussed. In Sec. V we take 
into account the cavity field damping and study its effect 
on motional state of the membrane. Finally, we summa- 
rize our conclusions in Sec. VI. 
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II. THE NONLINEAR COUPLING OF 
RADIATION PRESSURE FIELD TO A 
MECHANICAL DEGREE OF FREEDOM 



A. Model Hamiltonian 



As is shown in Fig.l, we consider a high finesse cav- 
ity which is detuned by the motion of a partially reflec- 
tive membrane placed between two macroscopic, rigid, 
perfectly reflecting fixed end mirrors [2TH231. [43] . Unlike 
the standard optomechanical structures, in this type of 
optomechanical system the coupling between the mid- 
dle membrane and the optical cavity strongly depends 
on the position of the membrane. This position de- 
pendence results in a cavity detuning, which is a pe- 
riodic function of the membrane displacement x, i.e., 
u! c (x) — (cj L)cos~ 1 [\r c \cos(Airx/ A)] where L and r c are 
the cavity length and the field reflectivity of the mem- 



brane, respectively. The Hamiltonian of the system is 
given by (excluding damping and driving terms) [21[ 122) 

H = hui c {x)a) a + HcOmb^b, (1) 

where a and b are, respectively, the lowering operators for 
the optical and mechanical modes satisfying the commu- 
tation relations [a, a^} — 1, [b, tf] — 1, and ui m is the oscil- 
lation frequency of the middle membrane. The fact that 
the motion of the membrane is quantized allows writing 
the operator of membrane's displacement x in terms of 

the Lamb-Dicke parameter 77 = ^/ 2 „!L iH US] (wp 
and to denote, respectively, the frequency of the incident 
field and the motional mass of the membrane)as 

c 

By expanding oj c (x) in term of b and w we obtain 



lj c (x) = 



7Tr 

2L 



2L 



E E 

m=l fc=0 



m 
k 



f r in6(m-2k)(b+tf) .1,-1 

-Iff ro-l yi2 1° T ft.o.j, 



(3) 



where to = 21 + 1 (I is an integer number) and 9 = 
2Luj m /c. By using the Baker-Campbell-Hausdorff theo- 
rem in Eq.(3| we may rewrite the Hamiltonian of equa- 
tion ([I]) in the form 



where 



H — H a + H 3 mt , 



Hq = huioa'a + HcOmPb, 



(4) 



(5) 



with wq = 7rc/(2L) as the natural frequency of the cavity 



without middle membrane, describes the free Hamilto- 
nian of the quantized cavity field and the free motion of 
the mechanical degree of freedom, and 

Hit = ha^a[ X *f 3 {h b )y + Xjftyfi&b)}, (j = 0, 1, 2, ..) 

(6) 

where 



2L V A V 2 



(7) 



and the nonlincarity function fj(fib) is defined 
by(appendix A) 



m)= E E ^ ("I) 4 ^ ) J e^ W ' ( - tt) '} j ^£i[W a ('"-2t) a ], ( 8 ) 

m=l fe=0 v / 2 



with fit, — w b and as the associated Laguerre polyno- 
mial, describes a nonlinear coupling of the radiation pres- 
sure field with the movable membrane through j-phonon 
excitations of the vibrational sideband. 



B. TIME EVOLUTION 

The time evolution operator corresponding to the 
Hamiltonian Q can be evaluated easily. For this pur- 
pose let us consider the first excitation of the vibrational 
sideband by choosing j = 1 in the Hamiltonian Q : 



H = hu! a^a + hui m b>b + h a' a[x* f {hb)b + xb fi^b)], 

(9) 
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where we have defined 

X = Xt 
f(n b ) 



2ttc 

i 

LX V 2mujr, 
fi(n b ). 



(10) 

The nonlinearity function fifib) plays an important role 
in our treatment since it determines the form of nonlin- 
earity of the intensity-dependent of the coupling between 
the cavity field and the membrane. We point out that in 
the limit of very small values of the Lamb-Dicke param- 
eter r\ and for certain values of the membrane reflectiv- 
ity r c the nonlinearity function f(rib) reduces to unity. 
Fig. 2(a) shows f(rib) as a function of rib for rj = 0.8 and 
r c = 0.99, while Fig. 2(b) displays the behavior of f(rib) 
versus n& for r\ = 10~ 5 and r c = 0.9. Obviously, for 
f(n>b) = 1 the Hamiltonian ^ reduces to the Hamilto- 
nian of the standard opto- mechanical system[57]. There- 
fore, the inherent nonlinearity of the model under con- 
sideration can be attributed to the parameters of rj and 
re- 
in the following, we consider the unitary time evolution 
operator in the interaction picture, corresponding to the 
Hamiltonian ^ that reads 



n J n int 



{t')dt'\ 



(11) 



where T indicates the time ordering operator. By sub- 
stituting Hamiltonian ([9]) in the interaction picture 



(12) 

into Eq.(|ll|) and using the Feynman operator calculus [46] 



to disentangling the time-ordered evolution opera- 
tor(appendix B) we obtain 



U. 



0=i) 



(r) 



e i\l3\ 2 \(T)(^a) 2 g(n b ) e O(l3\n b ) + ... 



(13) 



where by definition 

t = oj m t; (i = x/w m ; 

and 



H(r) = l-e* r ; X(T)=t + ifj*(T), 

(14) 



B = f(n b )b, 
g(n b ) = (fib + l)f 2 (n b ) 



(15) 



fibf 2 (n b - !)• 



Now we make use of an approximation to simplify the 
unitary operator in Eq.( |13[ ). In this approximation we 
keep terms up to second order in /3. The accuracy of our 
approximations may be verified by considering the real- 
istic values of j3 based on the experimental values of the 
relevant parameters [2 ll [2"2"]. For w m /27r = 10 5 Hz, m — 
50pg, X — 532nm and L = 0.67cm, we have \x\ ~ 10kHz 
which leads to ~ 0.01. Within this approximation the 
time evolution operator in Eq.(13) reduces to 



0=i) 



(r) 



e f3n{T)at a& e -t3*^(r)^aB e t\f3\ 2 \(T)(a^a) 2 g(f lb ) _ 



(16) 



Now, we assume that the membrane is initially pre- 
pared in a coherent superposition of the phononic Fock 
states \ip(0)) m — J2Ck\k} m and the cavity field is ini- 

k 

tially in a coherent superposition of the photonic Fock 
states \ip(0)) f = D n \n) t. Thus the initial state vector 

n 

of the system can be written as 



|^(0)) = EC^n|A;)m®|»>/. (17) 
k,n ^ ^ 

The time dependent state of the system is obtained by 



mr)) = ut 1] (r)\m) 



(18) 



By substituting (16) and (17) into Eq.(18) we obtain 



lV>(r)) - E C k D n e in2@ ^e A ^e- A nB\ k)m g, | n) /; 

n,k 

(19) 

where A„ = n^fi(r) and Q = |/3| 2 A(r). Thus, after some 
rearrangements we obtain 

IV-(t)) = E C k D n e™ 2& 9W\n) f ® |A n , fe (r)) m , (20) 



n , k 



where 



|A n>fe (r)) m = /(* - 1)!^: (A?l) I, A " y >< (^1 ' 



Ll' 



lll'l 



' (fc + QI f( k + l-l'-l)\ , 
(k-l')\ [/(fc-Z'-l)!] 2 1 + 



with k > I'. As is seen from Eqs. pOl and (21), if the 



membrane is initially prepared in a vibrational vacuum 
state, i.e., if C k = <5fc,o then after a time r = 27r it returns 
to its original state. At all times between r = and r = 
2-7T the membrane state is entangled with the cavity field 
state. Furthermore, we see from equation ( 20 ) that there 



is an explicit Kerr-like term in the time evolved state, so 
that physically one might expect the cavity field to have 
an evolution similar to that in a Kerr-like nonlinearity. 



III. NLCS OF THE MOVABLE MEMBRANE 
AND ITS NONCLASSICAL PROPERTIES 

A. Generation of NLCS of the membrane 

In this subsection, we consider the generation of NLCS 
of a macroscopic membrane. For this purpose we turn our 
attention to the state ( 20 ) . We assume a simple situation 
Cfe = <5fc,o that is the membrane is initially prepared in its 



motional ground state. Experimentally, the vacuum state 
of the membrane is feasible by cooling the membrane to 
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its ground state 
( 20 ) reduces to 



Thus, the time dependent state in which the state vector |Afc(r)) m (given by Eq.(23))is 

the right-hand eigenstate of the deformed operator C 
given by Eq.@ 



2 ^)\n) f ®\K{r)) m , (22 ) 



where |A„(r)) m = |A n0 (T)) m . The explicit expression 
for the state |A n (r)) m in the number representation is 
given by 



|An(r)>, 



„A„B 



t |0>m = ^|/) m m<«|A„(r)} m , (23) 



C\A k (r)) m = A fc (T)|A fc (T)), 



(31) 



In order to investigate the nonclassical behaviour of the 
generated NLCS, we focus our attention on the temporal 
evolution of the quadrature squeezing and of the Mandel 
parameter in the next subsection. 



where 



(I I An(r)>, 



K^fA V 



and the normalization constant reads as 



E 

i 



\pmi 



|A r , 



(24) 



(25) 



i-i 



with P(l) = f(l - 1)! = I] /(.?), (' > 0) and P(Q) - 1. 

3=0 

The density matrix of the system can be obtained by 
using Eq.@, 

p(r) = |V(t))(^(t)| 



B. Quadrature squeezing of the NLCS of the 
membrane 



Squeezed states are known by reduced quantum fluctu- 
ations in one quadrature of the field at the expense of the 
increased fluctuations in the other quadrature. In order 
to investigate the quadrature squeezing of the generated 
nonlinear coherent state of Eq. ( 23 ) we define 

1 



^(^-(Se^ + Ste-*-), 



£ 

n . I 



X 2 (r)=i(6 e --6t e --). (32) 

^ 26 These quadrature operators show the position and mo- 
D n D*e % ( n 6 ' ^ g ^\n)ff(l\ ® |A n (r)) mm (Aj(r)|m€aitum operators of the mirror, respectively. They obey 

the commutation relation 



and the reduced density matrix of the motional state of 
the membrane reads as 



Pm{r) = Trf[p(r)] 

= £p 



\K(r)) r 



(27) 
.<A„(r)|, 



where |A„(r)) m can be identified as a family of NLCSs 
of the membrane at time r whose amplitude depends on 
n. In order to get more clear insight to the above result 
we present the following argument. 

In principle, the nonlinear coherent states |C;/) are 
defined as the right hand eigenstates of the general- 
ized (deformed) annihilation operator B = f(n)b[. 



[X 1 (r),X 2 {T)]=i/2, 
and, consequently, the variances 
< (AA^-(t)) 2 >^< {X 3 {r)f >-< (X 3 (r)) >\ 
satisfy the uncertainty relation 



((AX 1 (t)) 2 )((AX 2 (t)) 2 )>1/16. 



(33) 



C? = l,2) 
(34) 



(35) 



A quantum state of the membrane is said to be squeezed 
when one of the quadrature components Xi and X 2 sat- 
isfies the relation 



£|C;/) = ClC;/>, 



(28) 



where f(n) is a reasonably well-behaved real function 
of the number operator h — wb and £ is an arbitrary 
complex number. We may seek for an operator & which 
is conjugate of the operator B, that is [B, C'] = 1 while 
their Hermitian conjugates W and C satisfy the dual 
algebra [C, B^} = 1. Thus it is easily found that 



C = 



' b, & = 5t ' 



(29) 



Now, let us assume D n = 5 n ^ in Eq.(22), i.e., the cavity 
field is initially prepared in \ip(0)) f = \k)f. Thus the 
reduced density matrix of the membrane reads as 

Pm (r) = e lfc2(e - e * )9(0) |A fc (r)) mm (A fc (r)|, (30) 



(Alj(r))^ < 1/4, (j = lor2). 



(36) 



The degree of squeezing can be measured by the squeez- 
ing parameter Sj(j = 1,2) defined by 



S J (r)E4<(Ai J (r)) 2 >-l. 



(37) 



Then, the condition for squeezing in the quadrature com- 
ponent can be simply written as Sj(r) < 0. The squeez- 
ing parameter Sj(J = 1,2) can be expressed in terms 
of the phonon annihilation and creation operators of the 
membrane as follows 



Si(r) = 2A 1 (t) + 2Re[A 2 {t)] - 4(i?e[A 3 (r)]) 2 , 
S 2 (t) = 2At(r) - 2Re[A 2 (r)} - 4(Im[A 3 (r)]) 2 , (38) 
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where 



IV. GENERATION OF SUPERPOSITION OF 
THE MEMBRANE'S NLCSs 



Ai(t) =< tfb >, A 2 (t) =< b 2 > e 2lT , A 3 (t) =<b> e". 

(39) 

Now, we study the temporal behaviour of S% (r), which 
gives information on quadrature squeezing of Xa(r). Nu- 
merical results are presented in Fig. 3, where we have plot- 
ted 52 (r) versus the scaled time r for different values 
of the parameters r c and 77. As is clear, the quadra- 
ture component X 2 (t) exhibits squeezing repeatedly at 
times r = (2m + l)ir, (m = 0,1,2,...) and the ampli- 
tude of squeezing is strongly influenced by the Lamb- 
Dicke parameter 77. Fig. 3(a) shows that for r c — 0.9 the 
quadrature squeezing only exists for 77 > 0.19 and the 
amplitude of squeezing increases by increasing the pa- 
rameter 77. Furthermore, we find that with the increasing 
value of membrane's reflectivity, the quadrature squeez- 
ing is strengthened and it appears for smaller values of 
ri. For example, Fig. 3(b) shows the quadrature squeez- 
ing of X 2 (r) for r c = 0.98(this relatively high value of the 
membrane's reflectivity has recently been considered in 
Ref.[47]). As is seen, for 77 > 0.14 the quadrature compo- 
nent X 2 {t) exhibits squeezing periodically in the course 
of time evolution. 



C. Mandel- Parameter 



The Mandel-Parameter corresponding to the NLCS of 
Eq. ( 23 1 is obtained as follows 



M(r) = - I 



E 



(Mr)) 

IT 

fe! 



(40) 



X 

k 



K 2 \A n {T)\ 2k kP{k) 2 \2 
fc! ) 



^ fe! 

k 



This parameter vanishes for Poisson distribution, is pos- 
itive for the super-Poisson distribution, and is negative 
for the sub-Poisson distribution. Fig. 4(a) shows the time 
evolution of the Mandel parameter with respect to the 
scaled time r for two values of 77 and for r c = 0.9. As 
is seen, the generated NLCS for the motional state of 
the membrane exhibits the sub-Poissonian statistics in 
the most of time. In addition, the sub-Poissonian char- 
acteristic is enhanced by increasing the Lamb-Dicke pa- 
rameter. Fig. 4(b) shows the behaviour of the Mandel 
parameter for a higher value of the membrane's reflectiv- 
ity r c = 0.98. As is clear, the higher values of r c leads to 
a considerable enhancement of the sub-Poissonian char- 
acteristic. 



A. Superposition of NLCSs with different 
amplitudes 

The system under consideration may also be used to 
generate a superposition of separated NLCSs of the mem- 
brane. In order to verify this claim, we assume that the 
cavity field is initially prepared in a coherent state, i.e., 

D n = e 



Thus the reduced density matrix of the 



membrane, Eq.(27), takes the following form 



2/i 



p m (r) = e -H- V ^e^(^ e *)^°)|A n (r)) mm (A„(r)|. 

n 

(41) 

This state shows a superposition of NLCSs with different 
amplitudes A„ for each value of n. The Q-function of the 



motional state associated with the state) 41 1, is given by 



Q(r) ee 



^<7|p m (r)| 7 ) = ie-l-l "M £g„, 



where 



(42) 



M 2 " in 2 eg(0)-m 2 6'n(0) 



E 

Z = l 



II 



(43) 

Fig. 5 illustrates the normalized Q-function(Q(T)) for dif- 
ferent values of the Lamb-Dicke parameter and the mem- 
brane's reflectivity at time r = 2.9 and for \a\ 2 — 4. As 
is seen from these figures the separation between nonlin- 
ear coherent components is increased by increasing the 
parameters 77 and r c . Figures. 5(a) and 5(b) show the 
double-peaked structure of the Q-function. We see that 
how the separation between the two components can be 
controlled by varying the parameters 77 and r c . In other 
words, by varying 77 and r c one can control the macro- 
scopic distinguishability of the states involved in the su- 
perposition. Fig. 5(c) shows an interesting case in which 
the membrane's state is a superposition of two different 
kind of states with equal amplitude, one of them is a co- 
herent like state and another is a squeezed like state. We 
mention that this state can only be prepared for large 
values of 77 = 0.98 and r c = 0.998. 



B. Nonlinear Schrodinger-like cat state 

In this subsection we consider the generation of the 
so-called nonlinear Schrodinger cat state (NLSCS) of the 
membrane. For this purpose, we assume that the mem- 
brane is initially prepared in the NLCS discussed in sec- 
tion. Ill, i.e., |r/>(0)> m = IC;/)m = ECfc|fc)m with 

fc 



C k = K 



/fc! 



, where the normalization constant is given 
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ICI 2fc -P(fc) 2 



Thus the time dependent state 



by W 

of the system given by Eq.(|20|) at time t = 2ir reads 



|^(T = 27r))=N'5^D n ( 



n. . k 



l)f |fe) T 



(44) 



It should be noted that the experimental realization of 
the system under consideration shows i] < 1 and 6 « 
l,(e,g. for the experimental values given in Refs. [2"TH2"3"] . 
we obtain 1CT 2 < r\ < 1 and 9 < 10~ 3 ). Therefore, one 
may keep terms up to first order in the phonon number 
fit, in Eq.Q and approximate the nonlinearity function 
f(j=i) i n b) = fi^b) by expanding the associated Laguerre 
polynomial 



where 



an b , 



(45) 



e= E E (m - 2fc) 

m=l fe=0 



m=l fc=0 



(46) 



Under this approximation, the nonlinear function g(rib) 
reduces to 



g(k) ~ Tfc 2 
where we have defined 



Afc 



r = (e - a) 2 + at + 2cr; A = a{a - 3e) + e 2 



(47) 



(48) 



By using Eq.(47) in Eq.(44), and assuming that the 



cavity field is initially prepared in the Fock state |/), 
D n = 5 nt i, the state of the membrane is obtained as 



\k)r 



where 



2nl 2 



(49) 



(50) 



Depending on the value of the parameter £, the state 
|'0(27r)) m can be made equivalent to a variety of nonlinear 
multicomponent cat states. As an example for, £ = 0.25 
we have 

| V(27r)> ro = e 2 "^\ + (e"/ 4 |C; f) m + e~ W4 | - C; f) m ) ■ 

(51) 

(-C 2 )"P(«) : 



- = L^2H' 2 E„ =0 



2 is the nor- 
The state fBTl) shows a two- 



where H+ = [2 + 2F' 2 
malization constant 
component NLSCS. The all higher order-multicomponent 
nonlinear cats states of the motional state of the mem- 
brane may be obtained at time r = 27r by varying the 
parameters £ and The normalized <5-functions(Q(T) = 
Qm(i~) of the cat states produced by different values 



of £ and £ and for a given value of the membrane's reflec- 
tivity are shown in Fig. 6. It is interesting to note that 
by varying the membrane's reflectivity one can control 
the separation between the components of the generated 
NLSCS. To show this, we have plotted the normalized 
Q- function for three different values of r c in Fig. 7. Ob- 
viously, by increasing the parameter r c , the separation 
between the components of the produced NLSCS is in- 
creased. 



V. EFFECTS OF CAVITY FIELD DAMPING ON 
THE MOTIONAL STATES OF THE MEMBRANE 

In order to make the model under consideration more 
realistic, we now take into account the cavity field damp- 
ing to examine its influence on the generated motional 
states of the membrane. For this purpose, we consider a 
regime in which the radiation mode relaxes much faster 
than the mirror. Experimentally, in the system under 
consideration the membrane Quality-factor is too large 
(Qm = I0 6 [H1 H2]) which leads to small membrane 
damping rate compare to the cavity damping k. The case 
in which the mirror relaxes much faster than the cavity 
mode, does not show any quantum features due to the 
thermalization effects [2S]. We assume that the number 
of thermal photons is negligible at optical frequencies. 
Hence, the master equation for the whole system reads 
as EE] 



i>(t) = -\p,H} + L{p), (52) 



with 



L(p) — n/2{2apa) — a) ap — pa) a), (53) 



where k denotes the rate of damping due to the pho- 
ton leakage from the cavity. In order to solve Eq.(52| 



we introduce the new density operator R a s p(r) = 
U(t)RU*(t). Here, U(t) is given by Eq.|l6|. By us- 
ing this definition we rewrite the Eq.(52|as 



R=UH{URU^)U^ = L{R). 



(54) 



We now assume the above differential equation has a so- 
lution as R — Rq + i?i , where Rq is a time independent 

operator, i.e.,i?o — 0(this operator corresponds to the 
free solution of Eq.(52|). While, Ri satisfies the follow- 



ing time dependent differential equation 



i?i = L(Rq + Ri). 



(55) 



Since the damping term L(p) is small enough, one can 
apply the first order Born approximation to re plac e 
L(R + Ri) by L(R ) in the right-hand side of Eq.(55|. 
Under this approximation, the following solution for the 
operator R is obtained 



(56) 



R(t) =R + L(R ;t)dt. 
Jo 

Hence, the density matrix p(t) is given by 
p(t) = /5 (r) +pi(r), 



(57) 



where 

p (r) = U(r)muHr), (58) 
is the undamped part of the density matrix given by 



Eq.(26) and the damping part pi(r) has the following 
form 

p 1 (r) = U(t)[[ L(R ;t)dt}UHr). (59) 
Jo 

To study the time evolution of the system in the presence 
of cavity damping we turn our attention to the density 
matrix of Eq. ( 26 ) in which the cavity field is assumed to 



be initially in a coherent superposition of Fock states, and 
the membrane is assumed to be initially prepared in the 



motional ground state. By using Eq.(53l the damping 
part pi(r) reads as 



Pi( T ) = K / [«( T - i )Po('r)a t (r - t)]dt 
Jo 



(60) 



-K/2[a t a / 5 (r) + po{r)a} a], 
where we have defined 

l{r-t) = U{t -t)all\T -t). (61) 



Substituting from Eq.( |B2[ ) into Eq.|6IJ) and applying the 
Baker-Campbell-Hausdorff formula together with using 
the Feynman disentangling theorem we arrive at 



~a(r -t) = ae' 3 ^ r - t )^ + e - ' 9 *' I * (r-t) ^e i| ' 9|2A(r ~ t)ff( ™ 6) . 

(62) 

By substituting po(r) from Eq.( 26 ) into Eqs.(60 1 we read- 
ily obtain 



p 1 (r) = nY J D n D* l e 1 ^ 



e-i 2 e')g{o) 
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(S(t - t)\n) f ® |A n (r)) m ) ® m (A i (r)|a t (r - t)) 



dt 



K/2[J2(n + l)D n Dte^ n e ~ l e ' )a ™\n) ff (l\ ® |A n , (r)) rom <A, j0 (r)|]. 



(63) 



r 



Finally by using Eq.(62) we find the following expression 
for pi (r) 



ol |/3| 2 A(r-t)g(fe) p - l |/3| 2 A*(T-t)3(fe') 



\n - l)ff(l - 1| ® |Ai,fc(T - t)) mm (A hk ,(T - t)\ 



- K /2[^(n + /)A l A*e l(n2e - z2e * )ff(0 V)//ai®|A^oW) mm (A ; , (r)|], 



dt 
(64) 



where by definition 

T k - k '(n-\ - r i(n 2 e~l 2 e*) (l (0) ^ 2 V^Vl(^(r)) k (^(r)) k ' P(k)Pjk') 

Vz v) ~ e 

(65) 



and |A nj fc(r)) m is given by Eq.(2I) 
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We are now in a position to examine the influence of 
the cavity held damping on the motional state of the 
membrane by using the density matrix of the Eq.(57). 



As before, we assume the membrane to be prepared in a 
vibrational vacuum state initially. In Fig. 8 we have plot- 
ted the time evolution of the Mandel parameter of the 
membrane motional state for various values of the rate of 
cavity field damping, the membrane's reflectivity and the 
Lamb- Dicke parameter. From Fig. 8(a) it is evident that 
unlike the undamped case (k = 0), at the initial stages of 
evolution the Mandel parameter increases which shows 
super-Poissonian statistics. However, as time goes on, 
this parameter decreases and the sub-Poissonian statis- 
tics occurs. The rate with which the enhancement of 
super-Poissonian (sub-Poissonian) statistics occurs is di- 
rectly proportional to the membrane's reflectivity r c ; the 
greater r c is, the more quickly the Mandel parameter in- 
creases (decreases). Furthermore, as in the undamped 
case, with the increasing value of r c the sub-Poissonian 
characteristic of the motional state of the membrane is 
strengthened. Also, Fig. 8(b) shows that by increasing the 
lamb-Dicke parameter rj the sub-Poissonian behaviour is 
enhanced. In Fig. 9 we have illustrated the time evolu- 
tion of the squeezing parameter S2 (t) in the presence of 
the cavity field damping. As is seen, the photon leak- 
age from the cavity suppresses the quadrature squeezing 
of the motional state of the membrane in the course of 
time evolution. However, it is evident that increasing 
the Lamb-Dicke parameter 77 (Fig. 9a) or the membrane's 
reflectivity r c (Fig. 9b) brings about noise reduction in 
the quadrature X2 and thus will increase the quadra- 
ture squeezing. Thus, we conclude that the destructive 
effect of the cavity field damping on the nonclassicality 
associated with the motional state of the membrane can 
effectively be controlled by changing the parameters r c 
and 77 . 

As another aspect of the influence of the cavity field 
damping on motional state of the membrane, we turn 
our attention to the NLSCS introduced in Sec. IV. For 



this purpose, we consider the density matrix of Eq.( 57 ) at 
time t = 2ir in which pq(2tt) corresponds to the nonlinear 



multicomponent state of Eq. ( 49 ) and p\ (2ir) is obtained 



from Eq.(57). As an example, the corresponding normal- 
ized quasiprobability Q-function(Q m (T) /tt 2 ) at r = 2ir 
for two different values of the cavity field damping rate 
k and for given values of the parameters r c , £, £ is illus- 
trated in Fig. 10. As is seen, for small value of K the Q- 
function consists of two well-separated peaks (Fig. 10(a)). 
By increasing the damping rate, one not only decreases 
the height of the two peaks but also decreases the spa- 
tial separation between the NLCS components involved 
in the superposition (Fig. 10(b)). In this connection, 
further numerical analysis reveals that by changing the 
parameters r c and rj one can control the spatial separa- 
tion between the components of the generated NLSCS 
(i.e., controlling the influence of decoherence) . It is also 
worth noting that the generated motional NLSCS pre- 
serves more or less its quantum coherence even in the 



presence of the cavity field damping. This survival of 
quantum coherence, which can be attributed to the non- 
linear character of the system under consideration, is in 
full agreement with the fact that deformed (nonlinear) 
coherent states superposition can be more robust against 
decoherence than the usual (nondeformed) Schrodinger 
cat states [35]. It has been shown [3J)] that deformed 
states, due to their nonlinear character, give rise to a 
more rich phase space structure, part of which can easier 
survive against decoherence. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have introduced a physical scheme 
that allows one to generate and control the nonclassical 
properties of motional nonlinear coherent states and their 
superpositions for an undamped vibrating micromechan- 
ical membrane inside an optical cavity. We have shown 
that if the cavity field is initially prepared in a Fock state 
the motional state of the membrane may evolve to a fam- 
ily of nonlinear coherent states. We have been interested 
in analyzing the nonclassical properties of the generated 
state of the membrane, including the quadrature squeez- 
ing and the sub-Poissonian statistics. In particular, we 
have found that the Lamb-Dicke parameter and the mem- 
brane's reflectivity lead to an enhancement of the non- 
classical properties. As we have seen, with increasing the 
Lamb-Dicke parameter and the membrane's reflectivity 
the sub-Poissonian behaviour and quadrature squeezing 
of the motional state of the membrane are considerably 
strengthend. In addition, the scheme offers the possibil- 
ity of generating various types of the so-called nonlinear 
multicomponent Schrodinger cat sates of the membrane. 
We have shown that the separation between nonlinear 
coherent components is increased by increasing the pa- 
rameters rj and r c . We have also extended our treatment 
to a more realistic situation in which the photon leakage 
from the cavity as a relevant source of decoherence is in- 
cluded and have examined its influence on the nonclassi- 
cal characteristics of the generated motional states of the 
membrane. We have shown that it is possible to control 
the effect of the cavity field damping on the nonclassical 
behaviour of the motional state of the membrane via the 
Lambe-Dicke parameter and the membrane's reflectivity. 
In particular, we have found that the generated motional 
NLSCSs of the membrane can be more robust against 
decoherence than the usual Schrodinger cat states. 

At the end, we would like to point out that our treat- 
ment is restricted to the case of one-phonon excitations 
(j=l) of the vibrational sideband. It is expected to gen- 
erate some other interesting nonclassical motional states 
of the membrane, e.g. nonlinear squeezed states, by con- 
sidering higher order of excitations. We hope to report 
on this issue in a forthcoming paper. 
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By applying the Baker-Campbell-Hausdorff theorem in 
Eq. ([3]) to disentangling the exponential terms and by us- 
ing the series expansion of each exponential term we ob- 
tain 



«c(2) = E ~ & £ E 



2L 2L £-1 L-i m \ h ] Am-l\( 
m=l k=0 V * ' LV 



^&n ^ m2{m - 2kf Y,{[m{m-m 2l+j §§vV +h.c.}. (Al) 



By making use of the relation 



Eq.( Al I can be written as 



^c(i) = f£ - Mir}6y{fi{n b )V + /i.e.}, (A3) 



(n fc -;)!' 



(A2) where 



f -tHA - V V Ir^T ( m \ (m-l)l r -^fae) 2 (m-2fc) 2 1 ft 6 !(m-2fc)j [ [- ( V 0) 2 (m-2k) 2 ] 1 (n b +j)< 

m=l fc=0 \ / 



(=0 



(A4) 



Rewriting the above equation in terms of the associated By using the Feynman operator calculus |46j we readily 
Laguerre polynomials, L*(«) = £ (^Q^|^t , leads obtain 

to Eq.§. r>0=1) 



} (t) = e^' a &B exp[-ig*ata J B'(s)e~ iuJmS ds], 



Appendix B: Disentangling of the time evolution 
operator by using Feynman operator calculus 



where 



(B4) 



The unitary time evolution operator for the case of By applying the Baker-Campbell-Hausdorff formula 
one-photon excitations (j = 1) can be written as 



B'(s)(s) = e -/3M(s)ataBt £ e/3Al ( s) ataBt _ ^ 

ipplying the Baker-Campbell-Hausd 
together with the commutation relation (B3|, Eq.(|B5|) 



(Bl) 



where 



U, 



u= 



1} (t) = fexp( - itfa J {g&e iu} ™ s + g* B e - iul ™ s )ds) , 



can be written as 

B'(s)(s) = B+ £ hl^er- 1 ^(-ly-^Ftn + i), 

m—1 i—0 

^ (B6) 

where k(s) — Pii(s)aJa and F(h) = nf 2 (h — 1). Finally, 
by using the Feynman disentangling theorem we arrive 



(B2) at 



with B = f(n)b as the nonlinear (deformed) annihilation 
operator obeying the following commutation relation 



£/-0'=l)( T ) _ e ^ t (r)ataB t e - / 3>*(r)a t aB x 
e i|/3| 2 A(r)(a+ a) 2 g(n b ) e O(J3 3 ,n b ) + ... _ 



(B7) 



[B,&] = (h+l)f(h)-nf(n-l). (B3) 
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FIGURE CAPTIONS component cat states of the membrane at time r = 2ir 

in the presence of cavity field damping for two different 
values of the damping rate K, and for r c — 0.95, 
Fig.l: A high-finesse optical cavity with two rigid £ = 1-8, C = 0.25: (a)« = 0.01w m , (b)« = 0.4w m . 
end mirrors and a dielectric membrane centered at the 
middle of cavity. 

Fig. 2: The nonlinearity function /(n&) func- 
tion of phonon number n.;, for: (a) r c = 0.99, rj = 0.8. 
(b) r c = 0.9, r) = 10~ 5 . Here we have set 8 = 10~ 4 . 



Fig. 3: Time evolution of the squeezing parameter 
^(t), corresponding to the generated NLCS (23) versus 
the scaled time r: (a) for r c = 0.9 and rj = 0.14 (black 
line), -q = 0.19 (blue line), and rj = 0.24(green line). 
(b)for r c = 0.98 and r\ = 0.1 (black line), r) = 0.14 
(blue line), and rj = 0.18 (green line). Here we have set 
L = 0.07m, m — 50pg. 



Fig. 4: Time evolution of the Mandel parameter 



M(r) corresponding to the NLCS (23) as a function of 
t: (a) r c = 0.9, (b) r c = 0.98 with r] = 0.25(black line), 
and rj = 0.3(red line). 

Fig. 5: The normalized Q-function(Q(r) = <5 m (r)/H 2 ) 



of the generated state (41) of the membrane at time 



t = 2.9: (a) r c = 0.95,7? = 0.8; (b) r c = 0.998, t] = 0.82; 
(c) r c = 0.998,r; = 0.98. 

Fig. 6: The normalized Q-function(<5(r) = Q m (r)/N?_) 
of the multi-component cat states of the membrane at 
time r = 2tt for different values of £ and £, and for 
r c = 0.95: (a)£ = 1.1, ( = 0.25; (b)$ = 1.8, C = 0.25; 

( C )e = i.8, c = -\/i- 

Fig. 7: The normalized Q-function(Q(r) = <5 m (r)/H^_) 
of the multi-component cat states of the membrane 
at time r = 2tt for different values of the membrane's 
reflectivity r c , and for £ = l/y/8, ( = 1.8: (a)r c = 0.8. 
(b)r c = 0.87. (c)r c = 0.99. 



Fig. 8: Time evolution of the Mandel parameter 
M(r) for different values of the cavity damping rate: 
(a) rj = 0.19, and r c = 0.93, k = 0(black line), 
r c = 0.93, K = lu) m (green line), r c — 0.95, /t = lw m (blue 
line); (b) r c = 0.95 and rj = 0.19, K = 0(black line), 
77 = 0.19, k = lcj m (blue line), rj — 0.16, k — lw m (red 
line). 

Fig. 9: Time evolution of the squeezing parame- 
ter S-2 (t) for different values of the cavity damping 
rate: (a) r c = 0.95 and rj = 0.19, n = 0(black line), 
rj = 0.19, k = lw m (blue line), rj — 0.16, k — lw m (red 
line); (b) rj = 0.16, and r c = 0.96, k — 0(black line), 
r c = 0.96, k = lw m (blue line), r c — 0.95, k — lw m (red 
line). 



Fig. 10: The normalized Q-function of the multi- 
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